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We calculate line shapes of correlation functions by use of 
complete diagonalization data of finite chains and analytical 
implications from conformal field theory, density of states, 
and Bethe ansatz. The numerical data have different finite 
size accuracy in case of the imaginary and real parts in the 
frequency and time representations of spin-correlation func- 
tions, respectively. The low temperature, conformally invari- 
ant regime crosses over at T* w 0.7J to a diffusive regime 
that in turn connects continuously to the high temperature, 
interacting fermion regime. The first moment sum rule is de- 
termined. 



I. INTRODUCTION 

Dynamical correlations characterize the spectral prop- 
erties of physical systems. They are accessible by a mul- 
titude of experimental setups. The access to dynamical 
correlation functions for physically relevant systems is 
usually difficult even in exactly solvable models.ua Dy- 
namical spin-correlation functions in Heisenberg chains 
have been ^widely studied numericallyoEl as well as 
analyticallyflia The comparison of numerical and ap- 
proximate analytical results for the purpose of ac- 
curacy conttoJ— has been used in various previous 
approaches .Il3~t£l 

Usually the focus lies on the imaginary part of the 
correlation functions. The real and the imaginary parts 
can be Kramers-Kronig transformed into each other and 
thus hold the same information. This is also true for 
the Fourier transform. The information that can be ex- 
tracted from finite systems accessible by exact diagonal- 
ization (ED) concerning the thermodynamic limit is lim- 
ited. The accuracy of the results is different for different 
representations. In the case of finite systems it proves 
thus useful to actually calculate all three representations 
to extrapolate to the thermodynamic limit .Ell 

The dynamical correlation functions became system 
size independent for high excitation energieaHI or, equiv- 
alently, on short time scales.ta While finite systems 
thus allow for the determination of correlation functions 
in the thermodynamic limit at high frequencies or on 
short time scales, field theoretical results describe their 
asymptotic behavior on long time scales or for small 



frequencies.Em The perspective of this paper is to com- 
bine the strongholds of both methods. 

The system to be discussed here is the one-dimensional 
antiferromagnetic Heisenberg model 

H = (J SfSf +1 + J SfSf +1 + J z S[Si +1 ) 
i 

+ J 2 ^ {SfSi +2 + SfSf +2 + S[Si +2 ) (1) 



with the superexchange integrals J and J2 between 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
magnetic ions, respectively, z-axis anisotropy J z and 
spin-1/2 operator components S" with v = x,y,z at 
site I. Energies will be given in units of the in plane 
exchange, i.e., J = 1. This Hamiltonian is relevant for 
the description of the magnetic systems, in, many quasi-, 
one-dimensional materials as Sr 2 Cu03 Cs2CuCi4,E2l 
KCuF 3 ,El or CuGe0 3 .El 

We focus on the spin-correlation function 



X(q,ito n ) = j- 



dr e— <5|(r)S*,(0)) 



(2) 



with Matsubara frequencies u n = 2im/ (3, inverse tem- 
perature (3 = 1/T, (&B = 1) Fourier transformed 
spin operators in interaction representation S*(t) — 
c -Ht e~ zql S[ e Hr , and number of sites L. In its ana- 



lytically continued form, where iu>„ 
it determines the structure factor 



uj - 



-ie with e — > 0, 



1 Imx(<7, u>) 



7T 1 



-0ui 



(3) 



relevant for neutron scattering experiments. 



A. Numerical methods 

For finite systems the correlation function can be cal- 
culated through the diagonalization of the spin Hamilto- 
nian in the spectral representation since cigenfunctions 
\n) and eigenvalues E n are known. All numerical results 
in this paper are obtained using periodic boundary con- 
ditions. Defining the matrix elements 



V nm = (n \Sg\m) 



(4) 



and the Boltzmann factor 



1 



fnrn{P) = ^ (e 



-/3E„ 



(5) 



where Z 
write 



Tr e ° H is the partition function, one can 



lim Wnm(/W 



(co + £„-£; m )2+ e 2 



Im\(q,u) = 7T y^/nm(/?) |Km| 2 <5(w + £71 ~ -En 



(6) 



(7) 



The corresponding real-time retarded spin-correlation 
function is obtained via a Fourier transformation as 



It can be shown that Eq. (||) gives very accurate re- 
sults for the real part of the correlation function if it is 
determined at the dominant spectral lines cDj i!3 



fnm(P)\V n 



E n — E m ) 



\E n - E m - w i | - Aw) (12) 



The regularization parameter Aw can be set to zero if 
only excitations at Qj are present (define 8(0) = 0). For 
Heisenberg chains at intermediate temperatures and fre- 
quencies a choice of Aoj = 0.1J yields reliable results. 
For higher frequencies the results for the real part of the 
correlation functions are free of finite size effects. 



x(q,t) = 



l 

27 



duj e 



^W^/nm(/5) \V n 



,2 e ^„ 



-E m )t 



(8) 



where 0(t) is the Heaviside function. 

To determine the correlation functions in frequency 
space we use the same methods as described in Ref. [L7] 
which we briefly summarize. At low temperatures small 
systems exhibit a small number of dominant spectral lines 
at frequencies Wj-wifch usually can be attributed to spe- 
cific excitations. The imaginary part of the correla- 
tion function is determined most accurately by "binning" 
the data as 



m.n 



fnm(fl) | Kim | 2 [0(uJ n7n — UJ ini 



) - 6(uj nm - t^ up )] 



~ sup ~ inf 



(9) 



For small systems at low temperatures the appropriate 
choice is such that the interval boundaries lie in the mid- 
dle between the dominant spectral lines: 



~ sup 



-inf 



: CO 



3+1 



(Ulj + U3 j+1 )/2. 



(10) 



If only the "dominant" spectral lines arc present and 
if those lines form a well defined continuum in the ther- 
modynamic limit, i e for L — » oo, Karbach, Miiller, and 
coworkers have showr£j'E£rE3 that Eq. (|?]) can be used, 
appropriately scaled to the thermodynamic limit, by in- 
troducing a density of states with respect to appropriate 
quantum numbers derived from Bethe ansatz. This leads 
to the following representation of the imaginary part of 
the correlation functionO 



Qi=E„ 



£ 



2* f nm (0)\V nm \' 



(11) 



The sum covers only values of n and m such that Uj — 
E n — E m . In Heisenberg chains this representation is only 
applicable at T = 0. 



B. Field theoretical preliminaries and 
transformations 

The correlations described by X (q,uj) Eq- @ are dom- 
inant at q — 7r reflecting the antiferromagnetic instability 
of the system. We will thus focus on this wave vector. 
For q 7r and J 2 — the spin-correlation function has 
been studied in detail with bosonization techniques by 
SchulztJ and hasJater been improved including logarith- 
mic correctionsO The result of conformal field theory 
for any two-point function with scaling dimension x in 
Euclidean space (r, r) at low temperature T iso 



XcftO. 1 ") = Xo 



7TT 



sinh 7rT ( - + it) sinh ttT I 



(13) 



where v denotes the velocity of the low lying spin excita- 
tions, and Xo ^ s some constant. The spin wave velocity 
for frustrated Heisenberg chains has been determined nu- 
merically as v = 0.5tt(1 - 1.12 J 2 ) for J 2 < 0.2411.E3 The 
Fourier representation in momentum q and frequency u> 
space with Imw > is 

XcFT (q, lo) = sm(nx) v l ~ 2x Xo (vrT) 2 ^ 2 (14) 
F ^-v(q-ir) ^ p {u} + v(q-n) 



with 



F x {k) 



2irT 



2ttT 



dX 



o (sinh X) x 



2 x " 1 r(l-x) 



T{x/2-ik/2) 
Y(l-x/2-ik/2)' 



(15) 



The value of the scaling dimension x depends on the 
strength of the interaction or the anisotropy in case of 
a spin chain. For the XY model we have x = 1, and for 
the isotropic Heisenberg chain x — 1/2. In the lower half 
plane Xcft(Sj w ) ^ s gi ven by Eq. ( |l4[ ) with to replaced by 



2 



From this representation we learn that the function on 
the right hand side of Eq. (|l4|) is analytic in a strip around 
the real axis with |Imw| < x2irT as long as T > 0; for 
T = we have 



0, for oj < v\q — 7r] 



v 2 (q 



else. 



(16) 



These analytical properties are shared by the structure 
factor Scft(<7> w ), which is related to ImxcFT^' w ) y i a 
Eq. (||), i.e., it is analytic in |Imu>| < x2itT for T > 0, 
and ScftU, w) ^ w 2 *- 2 for T = 0. 

Performing the Fourier transform to real time we see 
that both Xc FT (<7,£) and Scft(<7, decay exponentially 
at finite temperatures and algebraically for T = and 
q = n, 



X CFT (n,t) 



exp(-a;27rT£), for T > , 7 ^ 
^~ 2x forT-0. 10 



For momenta q ^ n the function Xcft(9> decays expo- 
nentially with time t for any T > as well as T = 0. 

There are additional contributions to x(q,w) and 
S^q, w) on the lattice that are singular at finite values 
of uj even for T > 0. These contributions have their ori- 
gin in the existence of the lattice which leads to finite 
energy bands with upper band edge singularities. There 
are no universal predictions like for the lower band edge 
governed by conformal field theory and described above. 
An exception, of course, is the XY spin model which can 
be mapped to free fermions. 

As we discuss in Sec. IC the case of the XY model 
suggests to assume that \(q,uj) is singular at a frequency 
A where the imaginary part diverges like 



lmx(q,to ± ie) 



±(A-oj) a , forcj<A, 



0, 



for oj > A. 



(18) 



The upper (lower) sign yields the retarded (advanced) 
correlation function. If not stated explicitly we discuss 
the retarded functions. The Kramers-Kronig transform 
yields the singularity of the real part 

, . / cot7ra(A - oj) a , forw<A, 

In the neighborhood of a — we have a logarithmic 
singularity^ 



Re X (q,u) = -ln|A- u\ 



(20) 



Regarding the time dependence we note that both func- 
tions x(Qi t) an d S(q, t) are dominated by the singularity 
at A and show long time asymptotics 



x (q,t) ~t- (1+a) exp(-iAt). 



(21) 



Since the operator S z is self adjoint Imx(7r, oj) is odd 
in oj. In general we thus set Imx(?,w) ~ sign(w) (A 2 — 
oj 2 ) a . The Fourier transform FT[x(?r, oj)] is consequently 



identical to twice the sine transform of Imx(7r, u>) and 
x(7r, t) is real. For the XY case with a = —1/2 we like 
to note more explicitly the qualitative result 



Imx(w ± ie) 



± 



sign(oj) 



for \w\ < A, 



. 0, for \oj\ > A. 

with Kramers-Kronig transform 



(22) 



9 arsinh-i/(^) 2 -l . . 

2 . , v > v { , for \oj\ < A, 



for M > A. 



2 arcsm-y/ — 



(23) 



For overcritical frustration0 J 2 > J c = 0.2411 the 
Heisenberg chain exhibits a gapped spectrum with a 
lower bound f2 g . Considering square root divergences 
at the lower and upper edge of the spectrum 

Imx(^ie) = ( ± 7P^i^T' f0rns<M<A ' 
I 0, else. 



(24) 



we obtain the Kramers-Kronig transform 



RexM = 



0. 

± 



sign (a;) 



for fi g < \oj\ < A, 
, else. 



(25) 



The upper band edge singularities and the resulting 
algebraic real-time asymptotics exist only at sufficiently 
low temperatures. At intermediate temperatures the up- 
per limit of the continuum yields an anti-symmetrizcd 
Lorentzian contribution. 



lmx(q,u) ^ L-{4>) ~ L + {4>) . 



where 



L± = 



r cos <f>— ( A ± oj) I sin < 



(26) 



(27) 



r 2 + (A± w ) 2 

Limiting < 4> < n/2 the real part is simply given by 

Re X (g, u) ~L^- n/2) + L+(0 - n/2) (28) 
and the Fourier transform reads 



x(Q,u) 



,-rt 



sin(At + 4>) 



(29) 



This temperature range will be referred to as "diffusive 
regime" . 

C. XY model 

We demonstrate the overlap of the accurate short time 
scale results from the exact diagonalization of finite sys- 
tems and the asymptotic behavior accessible by field the- 
ory for an exactly solvable case, the XY model, where 



3 
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FIG. 1. Susceptibility for the XY model: (a) imaginary 
part, (b) real part. The full line in (a) is the exact results from 
Eq. the symbols are obtained via Eq. (pl|), and the step 

functions via Eq. (^). The dashed line is the result from field 
theory with upper band edge cutoff but without divergence. 
The lines in (b) are the Kramers Kronig transforms of the 
imaginary part, and the symbols are obtained using Eq. 
with Alu — 0. 



J2 = J z = 0. The spin operators in this model can 
be transformed to non- interacting- spinless fermions via 
a Jordan- Wigner transformation.Ej The structure ijactor 
Eq. can be given for L — > 00 in closed formfl The 
imaginary part of the susceptibility at q = it is 



Imx xy (ir, uS) — tanh(/3w/4) (4 - cu 2 ) 



(30) 



This is the field theoretical result Eq. ("J]) with scaling 
dimension x = 1 multiplied with the square root diver- 
gence at the upper band edge. The limit of T — ► is 
given by Eqs. @ and (§|). 

Fig. [l](a) shows the imaginary part of the susceptibility. 
The full line represents the exact results, the symbols are 
obtained via the regularization Eq. (fllj), and the step 
functions are given by Eq. (^). The dashed line is the 
result from field theory with an UV-cutoff but without 
the upper band edge divergence. Fig. 0(b) shows the real 
part, the lines are the Kramers Kronig transforms of the 
imaginary part, and the symbols are obtained using Eq. 
@ with Auj = 0. 

We conclude that the multiplicative approach of the 
low energy description from field theory with the high 
energy behavior is adequate. Also, the numerical ap- 
proaches give a reasonable approximation to the exact 
result. The values of the real part for uj > A show only 
very little finite size effects. The divergences of the real 
and the imaginary part at the upper band edge show the 
correspondence predicted by Eqs. and (23). 

The retarded, real-time correlation function can be de- 
termined numerically in the thermodynamic limit. 




FIG. 2. Real-time spin-correlation function in the XY 
case at T — (a) and T = 0.3 (b). Thermodynamic limit 
results (full lines) are obtained via Eq. The different 

broken lines show the deviation of results for finite systems 
Eq. (^). L = 14 yields a good representation of the correlation 
function up to t « 2/ J at T = and up tot w 8/J at T = 0.3. 
The thin lines show the asymptotic behavior as predicted in 
Sec. " 



It 



i6(t) lim 

L — *oc 



Tk+q) 



pi{Ek — Ek+ q )t 



(31) 



The energy dispersion is given by = Jcosfc, are 
Fermi distribution functions, and the sum covers the first 
Brillouin zone. In general for L > 10 the result is inde- 
pendent of L for all practical purposes. 

In Fig. U we show the retarded spin-correlation func- 
tion for finite systems compared with the result in the 
thermodynamic limit (L — > 00, full lines) at T — (a) 
and at T — 0.3 (b). The different broken lines show the 
deviation of results for finite systems Eq. (|J). L = 14 
yields a good representation of the correlation function 
up to t w 2/J at T = and up to t w 8/J at T = 0.3. 

The thin solid lines in Fig. gshow the asymptotic <~ a5 
behavior from the upper band edge divergence. The con- 
tribution from the low frequencies yields an additive term 
- t- 1 at T = and - e -^Tt with ^ — 1 a t 7 1 = 0.3 
as predicted in Sec. [F^. The fits have been obtained 
for 80 < t J < 100 to assure the asymptotic limit. The 
numerical data for finite systems and for T « 0.3 yield 
a good representation of the correlation function in the 
thermodynamic limit up to time scales that are already 
dominated by the asymptotic, large time scale behavior. 



D. Technical outline of the approach 

Renormalization group studies show that the XY 
model is one point of the line of critical fixed points to- 
wards which the interaction flows in a bosonized repre- 



4 



sentation of the Heisenberg model (J z — 1).0@ One thus 
expects qualitatively similar results for the unfrustrated 
Heisenberg chain as in the XY model. This should also 
hold for frustrated Heisenbeca-Ghains, at least for under- 
critical J 2 < J c = 0.241 lBEaH 

The discussion of the XY model implies that the rep- 
resentation of the imaginary part of the correlation func- 
tion is best achieved by multiplying the upper band edge 
behavior to the field theoretical expression. 



Imx(7r,w) = Im[x CFT (7r,w)] 



(A 2 - u> 



2 A 



,,, tf(A-H) 

(32) 



The real part of the susceptibility is given by the nu- 
merical Kramers Kronig transform (KKT) of Eq. (32). 
The real-time representation is obtained by the Fourier 
transform (FT) of x(7T, u>). We attempt this approach 
for frustrated Heisenberg chains, where the exact form of 
the correlation function is not known. 

The field theoretical expression in Eq. ( |32"| ) depends on 
the parameter of the scaling dimension x and a global 
prefactor v l ~ 2x \q- The latter is determined by request- 
ing the sum rule of the first moment of the susceptibility 
to be correctly reproduced. The sum rule 



1 



dcu to Imx(q', u>) 



(33) 



can be extracted very accurately from the finite size data 
as discussed in Sec. O. We are then left with the param- 
eter vector p(T) = [x, A, a] which we find to be temper- 
ature dependent. The scaling variable x determines the 
low frequency behavior of the imaginary and the real part 
of the susceptibility (Eq. (|l6|)) as well as the decay in real 
time space as given by Eq. (17). The upper continuum 
edge A positions the cusp or divergence of the real part 
(Eq. (|19|)) as well as the oscillatory behavior of x(ir,t) 
as a function of time (Eq. ([H|)). Finally the exponent a 
describes the shape of the real part cusp or divergence 
(Eq. (|l9|)) and the decay of the oscillations in time (Eq. 

©)■ 



II. SUM RULES AND PREFACTORS 

Since the imaginary and real representation of the spin- 
correlation function are Kramers Kronig related and, as 
a consequence of the bounded excitation spectrumEZI it 
is straight forward to find that the sum rule Eq. (^) is 
given by 



h(q,T) = - lima 



Rex(g,a;). 



(34) 



For limy^o h(Qt T) = 4AT 1 (q) the structure factor sum 
rule discussed in Ref. [l^ is reproduced. For arbitrary 
frustration J 2 the structure factor sum rules are con- 
nected to the ground state energy Eq of the system via 
#i(tt)/2 + Kx(0.5w) = 2E G /3. We recall that for J 2 = 
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FIG. 3. First moment of the spin susceptibility as ex- 
tracted from the asymptotics of the real part in finite sys- 
tems, (a) low temperature Enite size effects and thermody- 
namic limit values as extracted from the Enite size scaling 
represented in (b). (c) Temperature dependence. 



and J 2 = 0.5 one has K \{q) — 2(1 — cos q)EqMi. For 
J 2 = 0.5 the average (Sf Sf +2 }T=o = vanishes .Ej 

The crucial point is that Re.x(g, u> — > oo) depends only 
weakly on the system size LjlJ Figure |(c) shows h(q,T) 
as a function of temperature for different system sizes and 
frustration parameters. In Fig. ||(a) it becomes obvious 
that the result for J 2 = and L — 14 is for all prac- 
tical purposes in the thermodynamic limit for T > 0.3. 
Analyzing Eq. ( |l3| ) one finds the correlation length to be 
£ = v/(2nxT). For x ~ 0.5 and v = 0.5vr(l - 1.12J 2 ) we 
find that the finite size effects are of the order of 10~ 3 
when the correlation length becomes £ ~ L. 

Figure |](b) shows the values of Ii(tt, 0) as a function 
of the system size. We determine the thermodynamic 
limit with the algebraic scaling function Ix(tt,0,L) — 
Ix(tt, 0, oo) + AoL~ v . Systems with Lmod4 = (open 
symbols in Fig. ||(b)) in general converge differently than 
systems with Lmod4 = 2 (full symbols). The two cases 
yield two values the difference of which serves as an error 
estimate. For J 2 = we find /i(tt,0,oo) = 1.1821(7) 
which is very close to the exact value of ii(7T, 0) = 
1.1817258.... For J 2 = 0.15 and J 2 = 0.2411 we 
find /i(tt,0) = 1.185(10) and h(ir,0) = 1.173(2), re- 
spectively. A value of Ii(tv, 0) = 1.1521(5) has been 
found for J 2 = 0.35. The result for J 2 = 0.5 with 
Ji (ir, 0) = 1.000(3) is extremely close to the exact value of 
1. In Fig. ||(a) the peculiar finite size effects for J 2 = 0.45 
become apparent. Obviously they result from the gap 
valuoH-H of fig w 0.12 being just in the temperature range 
where the finite size effects appear. The determination 
of Ji(tt, 0) = 1.075(10) is thus less accurate. The results 
are represented by the symbols in Fig. ||(a). 
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In Fig. |3j(c) the temperature dependence of I\ (tt, T) is 
shown for different values of the frustration parameter. 
The general asymptotic behavior of lim-r-^oo xC 71 ") ^) ~ 
T -1 becomes apparent from the discussion in section 
HIE. For the first moment we find limj'_ ) . 0o I\ (n, T) 
0.5 /*, This is reminiscent of the structure factor sum 
rule 



lim Tlo Re y(tt, lo) = lim 



do' a/ 2 SfacS) = 0.5 



(35) 

and is generic for all values of Ji and the model. 
Note that lim/r^oo S(q, lo) — UruT— ►(» S(q, —u>). 

The values of Re%(7r, lo = 0) show little finite size ef- 
fects even at rather low temperatures. Figure |i] shows 
a Log-Log plot from Eq. (|^) for different values of the 
frustration J2 and chain lengths as a function of temper- 
ature (broken lines). The full line shows the asymptotic 
behavior 

/DO 
dJ S(n,u') = 0.25 
-00 

(36) 

which reproduces a structure factor sum rule and is iden- 
tical for all values of the frustration and the XY model. 
The field-theoretical prediction Eq. ( |l4| ) for the prefactor 
j 1 " 2 *- 2 w ith constant scaling variable x is clearly inappro- 
priate for the temperature range shown. Our analysis in 
Sec. Ill shows that the deviation results from an explicit 
temperature dependence of the x(T) as well as from the 
temperature dependence of the singularity at the upper 
band edge, i.e., A(T) and a(T). 

The inset of Fig. ^ shows the finite size effects at low 
temperatures. liniT^o Rex^: oj = 0) diverges for J 2 < 
0.2411 and saturates for J% > 0.2411 which is reminiscent 
of the presence of a gapO 



the temperature de- 
) for different values 



For completeness we show in Fig. f 
pendence of Rex(7r, 3.8) from Eq. (£_ 
of J2. The finite size effects are < 0.1% and hardly visi- 
ble on this scale (full lines L — 14, broken lines L = 12). 
We do not show plots for J 2 — 0.5 since the presence 
of bound states makes the result unreliable, c.f. section 



[IIC. 



III. FRUSTRATED HEISENBERG CHAINS 

We now turn to the determination of line shapes of the 
spin-correlation function in frustrated Heisenberg chains 
making use of the precise results obtained above. 

A. Critical frustration 

We first discuss the values of J z — 1 and J2 = J c at the 
quantum critical point making the field-theoretical re- 
sults eligible for comparison. In frequency space at T = 




1.0 

Tin 

FIG. 4. Log-Log plot ofRex(TT,iu = 0) from Eq. (§) as a 
function of temperature. The full line is the universal large T 
asymptotic result ~ 0.25/T. The inset shows the finite size 
effects at low temperatures, limr^o R.ex( 7r , w = 0) diverges 
for J 2 < 0.2411 and saturates for J 2 > 0.2411. 
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FIG. 5. Temperature dependence of Re x( 7r > 3.8) from Eq. 
(^). The finite size effects are < 0.1% and hardly visible on 
this scale (full lines L — 14, broken lines L = 12). 



for a 14 site chain there are four spectral lines at frequen- 



e W 



(14) 



{0.264, 1.309, 2.112, 2.437} which, by 



analogy to the dimer-dimer correlation functions ,£3 may 
be identified as the triplet excitations out of the ground 
state.cB It is thus reasonable to suppose them to form a 
well defined continuum in the thermodynamic limit and 
thus Eqs. ( pT| ) and (|l|) can be applied with Alo = 0. 

The imaginary part of the spin-correlation function is 
shown in Fig. ^|(a). The bins are obtained using Eq. (||) 
and the symbols are given using Eq. ([ll]). The symbols 
in Fig. ||(b) show the real part as given by Eq. (p^). 
For lo = finite size effects are significant since one ex- 
pects from field theory and by analogy to the XY model 
Rex^ T=0 ^(7r, lo = 0) — * 00. For lo > 3 the numerical re- 
sults are essentially in the thermodynamic limit as can 
be seen from the different broken lines in the inset of Fig. 
H(b) . The real-time representation of the spin-correlation 
function at T = is given in Fig. 0(a). As in the XY case 
the system with L = 14 yields a useful representation of 
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m[J] 

FIG. 6. Imaginary part (a) and real part (b) of the 
spin-correlation function in the frustrated Heisenberg chain 
at T — with J2 = 0.2411. The imaginary part for Enite 
systems is binned (Eq. M )), each bin holds one spectral line, 
symbols are from Eq. (UjD. The symbols for the real part are 
obtained by using Eq.(U2). The full lines are the theoreti- 
cal result from Eq. ( pi} ) and its KKT. Inset: enlargement of 
cutoff region with Suite size results from Eq. 
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FIG. 7. Real-time spin-correlation function for 

J 2 = 0.2411 at (a) T = and (b) T = 0.3. Broken lines 
are finite size data from Eq. (19), full lines are FT of Eq. 



the correlation function up to t 2/ J. 

The fit with the theoretical predictions from Eq. (p2j), 
its KKT, and FT are given by the full lines in FigsTg 
(a), 0(b), and[7|(a) using the parameter set po.24ii(0) = 
[0.50(1), 2.6(1), -0.10(7)]. The cutoff 2.5 < A is bound 
by the highest spectral lines-, which must lie in the 
continuum. Previous resultsGI suggest A to decrease 
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FIG. 8. Imaginary part (a) and real part (b) of the 
spin-correlation function in the frustrated Heisenberg chain 
at T = 0.3 with J2 = 0.2411. The imaginary part for hnite 
systems is binned (Eq. (0)^), the symbols for the real part are 
obtained using Eq. (HS). The full lines are the theoretical re- 
sult from Eq. ( |3^ ) and its KKT. Inset: enlargement of cutoff 
region with hnite size results from Eq. m). 



monotonously with temperature and for reasons of con- 
sistency A < 2.7. The three parameters are then de- 
termined by matching the first maximum as well as the 
slope for 1 < t J < 2 in the real time representation and 
the value of the real part for lo ~ 3.8 (c.f. Fig. (||)). The 
finite size effects require to allow for rather large error 
margins. 

The result of x is consistent with the prediction from 
field theory. The value of a ^ suggests a more com- 
plicated upper continuum edge than a simple ultraviolet 
cutoff. We emphasize that the overall prefactor of the fit 
function is fixed by the sum rule Eq. (|33j) and that values 
for A and a have been obtained without using the not so 
well defined binned data of the imaginary part. 

The plot of the real-time representation of the spin- 
correlation function at T = 0.3 in Fig. ^(b) reveals 
the temperature dependence of the parameter vector 
Po.24ii(0.3) = [0.48(1), 2.5(1), -0.25(5)]. The result for 
L = 14 yields a useful representation of the correlation 
function up to t ~ 4/ J. The full line is the fit from 
the FT of Eq. (|3^). The exponential fall off predicted in 
Sec. IB is confirmed and renders the value of the scaling 
dimension. The oscillations are much less damped than 
would be obtained with an upper band edge exponent of 
a = thus yielding the negative value of a = —0.25. The 
cutoff A is given via the period of the oscillations. 

Fig. ||(a) and ||(b) show the imaginary and real part 
of the frequency representation of the spin-correlation 
function at T = 0.3, respectively. The binned data 
for the imaginary part via Eq. (^) and the symbols 
for the real part via Eq. ([lj) are obtained for the 
sets of dominant spectral lines, in the case of L = 



7 



14 they are given by the frequencies Qj £ VV2411 = 
{0.311, 0.971, 1.517, 2.068, 2.301}. The condition of a well 
defined continuum with respect to Bethe ansatz quantum 
numbers is violated and thus Eq. ( pd| ) cannot be applied 
any more. For the real part the data are regularized with 
Ao; = 0.1, which is determined to give reliable results 
analogously to the dimer-dimer correlation functions.EZI 
An exception is made at ui = 0, where no regularization 
is applied (Auj = 0). 

The good correspondence of the field-theoretical fits 
from Eq. © and its KKT (solid lines) in Figs. | (a) 
and ||(b) proves the reliability of the parameters ex- 
tracted from the real-time representation. Especially 
the good agreement of the values of Rexi^, 0) and of 
Rex(7r,w > 2.7) (inset Fig. ||(b)) are non-trivial con- 
sistency checks. We expect a thermal smearing out of 
the small divergence at the upper band edge of which 
the shape is not known and which we did not account for 
(solid line in Fig. ||(a)). This might have a small influence 
on the parameters extracted and thus we adapted rather 
conservative error bars. The finite size data in both the 
real and imaginary part suggest a steeper slope in the low 
frequency dependence of the fitted curves. Together with 
the temperature dependence of the scaling dimension x 
this indicates the breakdown of the scale invaciance pre- 
dicted by field theory at finite temperatures .E3 



B. Unfrustrated Heisenberg chain 

Heisenberg chains without frustration are relevant for 
most of the magnetically quasi one-dimensional systems 
studied experimentally. Since the system is integrable the 
numerical data can be compared to results from Bethe 
ansatz. 

At T = the four spectral lines of the triplet excita- 
tions out of the ground state for a 14 site chain are at 



(14) 



= {0.307, 1.57, 2.56, 3.10)110 It 



frequencies ujj 6 W, 
is thus reasonable to suppose them to form a well de- 
fined continuum in the thermodynamic limit and thus 
Eqs. ( pd| ) and (12) can be applied with Alu = 0. Binned 
data for the imaginary part are obtained via Eq. (||). 

The imaginary part of the spin-correlation function is 
shown in Fig. ||(a). The real part in Fig. || (b) shows for 
lu = significant finite size effects since Rex^ T=0 ^ (n, lu — 
0) — > 00. For lu > 3.5 the numerical results are essentially 
in the thermodynamic limit (Inset of Fig. |^(b)). The 
real-time representation of the spin-correlation function 
at T = is given in Figure |ic](a) . 

The upper edge of— the two-spinon continuum is know 
exactly to be A = nE3 Bethe ansatz results suggest that 
the infra red divergence of the two-spinon contribution 
%( 2 ) to the imaginary part of the spin-correlation function 
has a logarithmic correction 



lim H ^ Imx (2) (7r,cj) cx lu 1 y/ln(u> *) 



(37) 




FIG. 9. Imaginary part (a) and real part (b) of the 
spin-correlation function in the unfrustrated Heisenberg chain 
at T = 0. The binned curves in (a) are from (Eq. (ISpj, the 
symbols are from Eq. pj|). The symbols for the real part (b) 
are obtained using Eq. (H3). The full lines are the fits from 
Eq. ( p^ ) and its KKT. Inset: enlargement of cutoff region 
with finite size results from Eq. (la). 




while at .-the upper continuum edge it vanishes square 
root likeJlj The two-spinon contribution has been found 



FIG. 10. Real-time spin-correlation function for J2 — at 
(a) T = and (b) T = 0.3. Broken lines are finite size data 
from Eq. (g), Ml lines are FT of Eq. 

to contribute 72.89% to the total spectral weight .111 Our 
and previous numerical studies show that the spectral 
weight of the total correlation function for above the 
two-snjaon continuum (lu > n) at q — -k is less than 
0.1%.qB Thus the total spin-correlation function includes 
also about 27% higher order contributions and we have 
Imx(7r,ti;) > Imx' 2 '(vr, ui). Consequently we must re- 
quire x < 0.5, a < 0.5, and A = w. 

From the amplitude in the real-time representation in 
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FIG. 11. Imaginary part (a) and real part (b) of the 
spin-correlation function in the unfrustrated Heisenberg chain 
at T — 0.3. The step function in (a) is from (Eq. (19) ), the 
symbols in (b) from Eq. @. The full lines are from Eq. @ 
and its KKT. Inset: enlargement of cutoff region with results 
from Eq. @. 
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FIG. 12. Imaginary part (a) and real part (b) of the 
spin-correlation function for J2 = 0.5 at T = 0. The binned 
curves in (a) are from (Eq. (0)J, the symbols are from Eq. 
(fn[). The symbols for the real part (b) are obtained using 
Eq. (|l|). The full lines in are the fits from Eq. @ and its 
KKT. 



Fig. (a) we find that the parameters x and a fall 
on a line defined by (0.38,0.25) < (x,a) < (0.44,0.5). 
Taking also the value of the real part at ui = 4 in 
the inset of Fig. |(b) into consideration we determine 
p (0) = [0.40(3), 7T ± 0.01,0.33(5)]. The error margins 
have been chosen rather large because of the obvious fi- 
nite size effects. The resulting fits with the theoretical 



predictions from Eq. (|3_2 
the full lines Figs. |9|(a), 



, its KKT, and FT are given by 
(b), and [l(j (a) and show satis- 
factory agreement with the results from finite systems. 

The short-dashed line in Fig. |l0|(b) shows the corre- 
lation function for 14 sites in the real-time represen- 
tation. The exponential decay in time imply a scal- 
ing dimension of x > 0.4 and the strongly damped 
oscillations an exponent a > 0.5. The real part at 
higher frequencies as shown in the inset of Figs. |TT|(b) 
requires the exponent to be a < 0.7. Together with 
a matching value for Rex^,^ — 0) the parameter set 



Po(0.3) 



0.45(2), tt± 0.1,0.64(10)] yields the best fits 



from Eq. (52|)as shown by the full lines in Figs. fll](a), 
|Tl|(b), and 10| (b). The overall agreement is satisfactory, 
for small frequencies the numerical data suggest a slightly 
altered functional dependence on the -frequency than is 
reproduced by the field-theoretical fit.E3 



C. Overcritical frustration 



In the case of overcritical frustration the spectrum of 
the spin chains acquires a gap ri ff B3'tiJ We discuss here 



the value of J 2 = 0. 
suits from literature.! 



«fr-better comparability with re- 
For a two particle (spinon) 



continuum one expects a density of states that diverges 
square root like both at the Lower as well as a^-,the up- 
per edge. Previous numericaB and variationals results 
suggest sharp maxima in the density of states just above 
the lower edge il g and just below the upper edge A of 
the continuum accompanied by a square root like vanish- 
ing at both edges. This can be understood in connection 
with bound-states being present close to the edge of the 
continuum.t3 

At T — there are more spectral lines present than in 
the case of critical and undercritical frustration. Some 
of therruare signatures of the bound states present in the 
system.E3 We still apply Eqs. ( |TT| ) and (||) to extract the 
imaginary part of the spin-correlation function as shown 
in Fig. |l2](a). The fluctuations of the results at the upper 
band edge are reminiscent of the fact that the bound 
states do not form a continuum in the thermodynamic 
limit. We do not attempt to refine the plot by extracting 
the bound state contributions since results of the real 
part and especially the real-time representation are more 
reliable anyway. 

Fig. |l2|(b) shows the real part from Eq. ( |l2"| ) with Alu = 
0.1 for 1.5 < u) < 2.3 and Alu = 0.001 else. The real-time 
representation of the spin-correlation function at T = 
is given in Fig. |l3|(a). 

In order to adapt Eq. ( |32| ) to the gapped spectrum 
we have to include the value of the gap energy into the 
functions F x in the field-theoretical expression Eq. (|l4| ) 
as 
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7 2 =0.5, 7=0 
x=0.45, a=-0.6 




FIG. 13. Real-time spin-correlation function for J2 = 0.5 
at (a) T = and (b) T = 0.3. Broken lines are finite size data 
from Eq. (0), Ml lines are FT of Eq. Eh. 



F, 



uj ± v(q — 7r) 



F. r 



' u±^n 2 g + v 2 {q-Tr) 2 



2irT 



(38) 



At T = this yields the correct divergence at the lower 
band edge Imx(7r, w) ~ [lu 2 — fl 2 ] x , c.f. Eq. (|l6|), while 
for finite temperatures Imx(7r, u>) ~ w. 

The full line in Fig. |l3|(a) shows the lit from the FT 
of Eq. (H) with a: = 0.45(3), 0, = 0.25,B A = 2.18(10), 
and a = —0.6(1). The slow damping of the oscillations 
requires the large value of \a\, there is a small frequency 
modulation stemming from the lower boundary Q g . The 
agreement with the finite size data (broken lines) for the 
finite size effect free time domain t J < 5 is not nearly as 
good as for the cases of lower frustration. We conclude 
that the simple functional form of Eq. (32) is insufficient. 
The square root dependence reported in literature^ must 
be included. A perturbative examination for small J z and 
J2 showsE3 the strong interplay between the spinon con- 
tinuum and bound states. Since a theoretical expression 
for the correct spectral weight distribution is not known 
we limit ourselves here to the observation that the small 
damping of the oscillations in Fig. |l3](a) requires a rather 
sharp increase in Imx(7T, w — ► A). 

The full lines in Fig. |Ej(a) and [0](b) show the fit of 
Eq. ( |32| ) and its KKT to the imaginary-part and real-part 
data, respectively. They are similar to the forms given by 
Eqs. (|24j) and ( p5| ) where x — 0.5 and a = —0.5. Here the 
parameters x and a have been adapted to roughly match 
the numerical values of Rex(?r, 0) and Re x(tt, 3.8). Es- 
pecially the discrepancies of the fits at lu ~ A indicate 
that a more involved fit function is necessary. 

Fig. |l3|(b) shows that also at T = 0.3 an accurate fit 
(full line) with the functional form of Eq. d3j|) to the finite 
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FIG. 14. Imaginary part (a) and real part (b) of the 
spin-correlation function for J2 = 0.5 at T = 0.3. The step 
function in (a) is from (Eq. ^)), the symbols in (b) from Eq. 
@. The full lines are from Eq. @ and its KKT. 



size data (broken lines) is not possible. The parameters 
for the approximate fit are x = 0.48, Q g — 0.25,Ej A = 
2.2, and a — —0.55. Correspondingly, the fits to the 
imaginary- and the real-part representations in Fig. |l4|(a) 
and |l4|(b) show inconsistencies with the numerical data. 



D. Intermediate temperatures 

At intermediate temperatures the interaction in the 
system is expected to broaden out all sharp features in 
the correlation functions. The onset of this effect is al- 
ready observed at T = 0.3 as discussed in the previous 
sections. At T = 0.7 the exponent of the upper con- 
tinuum edge for J 2 = is a = 3.2(1) so that the sin- 
gularity is basically completely damped out. Also, the 
rather large effective continuum edge A = 3.45(10) does 
not quite reproduce the correct oscillatory behavior as a 
function of time. 

At about the same temperature the scaling dimension 
increases to x ~ 1. T* w 0.7 thus marks the crossover 
temperature from strongly interacting, conformally in- 
variant to noninteracting fermion and high energy diffu- 
sive behavior. This is consistent with c^Rex^, 00) and 
9-rRe y(7r, 3.8) being extremal at T T* as seen in Figs. 
I and § 

Figure |l5|(c) shows the real-time representation of the 
spin-correlation function of the unfrustrated Heisenberg 
chain at T — 1. The amplitude of the modulations be- 
tween tJ~2 and the onset of finite size effect for the 14 
site chain at t J ~ 7 cannot be fitted algebraically. The 
exponential fit from Eq. (E9h shown by the full line in 



Fig. 15(c) matches excellently. The analogy to the XY 



model suggests that the long time asymptotics is cap- 
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FIG. 15. Spin-correlation function of the unfrustrated 
Heisenberg chain at T = 1. The broken line in (a) is the 
binned imaginary representation, the broken lines in (b) are 
the real part from Eq (Q) with two values of e for regular- 
ization, all for L = 14. Full lines in (a) and (b) are double 
Lorentzian fits, (c) shows the real-time representation from 
Enite systems (broken lines) and the asymptotic fit (full lino) 
with A = 2.31(1), T = 0.941(5), and A = -0.67(5). 



tured in this fit. The extracted value for A = 2.31(1) 
thus marks the effective, thermally smeared out upper 
band edge. For the appropriate interpretation of the pa- 
rameter T = 0.941(5) refer to section IV. 

The discrepancy of the real-time fit function (full line 
in Fig. |l5](c)) and the correct line shape for small times 
does not allow for a direct comparison of the results 
with its Fourier transforms. The difference between the 
fit and the exact result is roughly exponential. The fit 
with Eq. (|2^) implies that the real and imaginary part 
should contain contributions from the continuum bound- 
ary Lorentzians Eqs. (^8|) and (p6|). Indeed, the double 
Lorentzian fits (full lines in Fig. |l5[(a) and |l5|(b)) with an 
additive Lorentzian contribution centered at ui = com- 
pare well with the binned data for the imaginary part 
and the real part data from Eq. (^) with e = 0.03 (bro- 
ken lines). The fit parameters even though similar are 
not such that the fits are appropriately Kramers-Kronig 
and Fourier related. The fits must thus be regarded as 
sophisticated guides to the eye. Similar results are ob- 
tained for J 2 > 0. Similar line shapes are also found in 
systems with large spins O 



E. High temperature limit 

In Fig. |l6|(a), |l6](b), and |l^(c) we show the respective 
imaginary, real, and real-time representation of the sus- 
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FIG. 16. Imaginary (a), real (b), and real-time (c) rep- 
resentation of the susceptibility for T — > oo. Broken lines in 
(a) are binned, in Eq. m) with e = 0.02 was applied, and in 
(c) Eq. (^). The solid line in (c) is the sinusoidal fit for the 
example of J 2 — 0. 



ceptibility for different values of the frustration in the 
limit of infinite temperatures. The data of the imaginary 
part (a) are binned and the real part (b) is given by Eq. 
(||) with e = 0.02 for L = 14. The time representations 
(c) from Eq. d§|) show finite size effects for t > 6/ J. 

For u> — > the slopes of the imaginary part are similar 
to the exact result for the XY model Eq. (p^p|). A small 
frustration dependence becomes obvious when plotting, 
the structure factor instead of the correlation functionrrEl 

The oscillations in time shown in Fig. |l^(c) can be 
fitted very accurately for 1.5 < tJ < 6 with an ex- 
ponential decay via Eq. (p9]). The parameter sets 
[A,T,4>] are obtained as [2.15(1), 0.375(2), -0.87(1)] for 
J 2 = 0, [2.31(1), 0.323(2), -1.13(1)] for J 2 = 0.2411, and 
[2.21(1), 0.345(2), -0.82(1)] for J 2 = 0.5. The full line 
shows the resulting fit function for J 2 = 0. 

In the classical limit, where (S^)t— «x> ~~ * °°, paramag- 
netic behavior is expected for T ^> J. This leads to an 
expected functional dependence of the structure factor 
of limx^oo Sdassiq, u) ~ lim e ^o e/(^ 2 + e 2 ). From Eq. 
@ follows that lim T _ K3O TRex(7T,0) = (5f) T ^oo which 
is consistent with the expected functional dependence in 
the classical limit. 

The XY model is one point of the line of critical fixed 
points towards which the interaction flows a bosonized 
representation of the Heisenberg modelOO The suscep- 
tibility of the XY model shows a square root divergence 
at us — 2 and Im^yfi^ > 2) = 0. The shape of the 
spectrum at the upper band edge observed for Heisen- 
berg chains is thus an interaction effect. □ The shape of 
Imx(7r,cj ~ A) indeed resembles that of a Fermi distri- 
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FIG. 17. Extracted values for (a) the sum rule Ii(ir, 0), (b) 
the scaling dimension x, (c) the upper edge of the continuum, 
and (d) the exponent at the upper edge of the continuum as 
a function of J2 at T = 0. 



bution of weakly interacting electrons. We thus interpret 
the limit T — > 00 as best described by weakly interacting 
spinless fermions. 



IV. RESULTS 

Figure [l7| summarizes as a function of the frustration 
J2 at T = the extracted values for (a) the sum rule 
h(n, 0), (b) the scaling dimension x, (c) the upper edge 
of the continuum A, and (d) the exponent at the upper 
edge of the continuum a. 

(a) The values of I\ for J2 < 0.2411 are within 
error bars almost identical which underlines the com- 
mon featu res , o f . H eisenberg chains with undercritical 
frustrationJiaO'Eafea 

(b) The scaling variable x shows a stronger infrared 
divergence for unfrustrated Heisenberg chains than for 
those with critical frustration. The values for overcrit- 
ical frustration have to b e regarded as effective ones as 
discussed in Section [II C . 

(c) The cutoff frequency of the upper limit of the 
spinon continuum is linear as a function of frustration 
for J 2 < 0.35. 

(d) The exponent of the cusp at the upper boundary of 
the spinon continuum is always smaller than the value of 
a = 0.5 predicted for the two-spinon contribution for J 2 . 
The value of a vanishespfor J 2 ~ 0.2 in agreement with 
the previous observations] that for that value the spectral 
properties of the frustrated Heisenberg chain are, similar 
to the conformally invariant Haldane-ShastryOc3 model. 

The prefactor xo from Eq. (|l4|) is of order 1 and 
slightly frustration dependent. The values (J2,Xo) 




FIG. 18. Extracted values for (a) the scaling dimension 
x for T < T* , (b) the upper continuum cutoff frequency for 
T <T* (full symbols) and the effective upper continuum edge 
in the diffusive/weakly-interacting- fermion regime (open sym- 
bols), (c) the cutoff exponent, and (d) the control parameter 
F determining the decay as a function of time in the diffu- 
sive/weakly-interacting- fermion regime. 



are: (0,1-31(5)), (0.15,1.12(5)), (0.2411,1.01(5)), and 
(0.35,0.88(5)). For J 2 = 0.5 one has v 1 ' 2 *^ = 0.69(5). 
Values for J 2 = 0.45 are not computed because of the 
peculiar finite size effects shown in Figure [| 

Figure [l8] summarizes the temperature dependence of 
the fit parameters for the experimentally most relevant 
unfrustrated chain with J2 = 0. 

(a) The scaling variable approaches the value of the 
XY model limit at the crossover temperature to the dif- 
fusive regime T* 0.7. The direct determination of 
x(T > T*) is not possible but since for T — > 00 the weakly 
interacting fermion case is recovered it is expected to lock 
in at x(T > T*) = 1. 

(b) The upper continuum edge A(T < T*) marks a 
sharp cutoff (full symbols) while in the diffusive regime 
(open symbols) it is the effective, thermally smeared out 
upper continuum boundary. In the weakly interacting 
fermion limit it saturates at A(T — > 00) = 2.15(1). 

(c) The exponent of the upper continuum edge a in- 
creases with increasing temperature reflecting the ther- 
mal smearing out of the singularity. Its value at T = T* 
is so large that the cutoff is barely singular and thus not 
very well defined. In the diffusive regime (T > T*) this 
quantity is undefined. 

(d) r is an effective parameter that controls the contin- 
uous transition of the system from the diffusive behavior 
at T* to the weakly interacting fermion limit at T — ► 00. 
For J 2 = it saturates at T(T ->oo) = 0.375(2) 
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V. CONCLUSIONS 

We have compared numerical results from the exact 
diagonalization of finite systems with results from con- 
formal field theory together with implications from the 
density of states, the exactly solvable XY model, and 
Bethe ansatz solutions for integrable systems. We use the 
different finite size accuracy of the imaginary, real, and 
time representations of the spin-correlation functions to 
extract reliable information on the thermodynamic limit. 

At low temperatures the dynamical correlation func- 
tions of frustrated Heisenberg chains are well described 
by a multiplicative superposition of the contribution from 
low lying elementary excitations described by conformal 
field theory and a density of states and matrix element 
induced singularity near the upper edge of the two-spinon 
continuum. At the frustration value of J2 ~ 0.2 the 
system is closest to the conformally invariant Haldane- 
Shastry model. 

At T* « 0.7 we observe the crossover from the low 
temperature, conformally invariant regime to a diffusive 
regime. All correlations in time then decay exponentially. 
The diffusive regime connects continuously to the weakly- 
interacting- fermion regime for T — > 00. 

We give the frustration dependence of the control pa- 
rameters for the line shapes of the spin-correlation func- 
tions at T — and their temperature dependence for the 
experimentally most relevant case of J2 = 0. The tem- 
perature dependence of the first moment sum rule of the 
spin-correlation function is accurately determined. 
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